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The use of novel perturbatively-improved gluonic actions on anisotropic lattices in which the temporal spacing 
is much smaller than that in the spatial directions is discussed. Such actions permit more efficient measurements 
of noisy correlation functions, such as glueball correlators, on coarse lattices. A derivation of these actions at 
tree-level is outlined; mean-field link renormalization plays a crucial role in their construction. Results for the 
low-lying glueball masses and the heavy-quark potential are presented. 



The use of improved actions makes possible ac- 
curate Monte Carlo simulations of QCD on coarse 
lattices with greatly reduced computational ef- 
fort. However, for some calculations, such as glue- 
ball masses 0, the coarseness of the temporal lat- 
tice spacing can be a severe drawback, greatly 
reducing the number of correlator time intervals 
which can be measured. This problem can be cir- 
cumvented by using anisotropic lattices in which 
the temporal spacing is much smaller than that in 
the spatial directions, enabling one to exploit the 
enhanced signal-to-noise of the correlation func- 
tions at smaller temporal separations. Improved 
lattice actions often admit spurious states with 
energies of order 1/a, where a is the lattice spac- 
ing. On fine grids, such modes are comfortably 
high-lying and little affect the low-lying modes 
of interest. However, on a coarse lattice, this 
may no longer be true. Reducing the temporal 
spacing also cures this problem by elevating the 
spurious states to higher energies of order 1/at, 
where at is the temporal lattice spacing. The lift- 
ing of spurious states is particularly important 
when designing improved fermion actions |^. In 
fact, anisotropic lattices will be useful whenever 
one is faced with a four-momentum in which one 
component is unusually large, such as in the cal- 
culation of glueball masses]^, the heavy quarko- 
nium spectrum, or hadronic form factors at large 
momentum transfers. Anisotropic lattices have 
also long been used in finite temperature studies. 

To remove O(a^) errors from all Green's func- 
tions in lattice QCD at tree level in perturbation 
theory, one needs only to eliminate O(a^) arti- 



facts in the lattice action. This can be done by 
performing a small-a expansion of some suitable 
lattice action and adjusting the interaction cou- 
plings so that the leading terms in the expansion 
reproduce the correct continuum QCD action and 
the 0{a?) terms are absent. Any lattice operator 
W^[C/] constructed from link variables will have 
a small-a expansion of the form 
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W,[U] = Y1 d'^C'a^'Y. ^i'^Qi'H^)' (1) 
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where Ca"^ are the expansion coefficients, the spa- 
tial lattice spacing is Og, the temporal spacing 
is at — ^Cs, and Qi^\x) are local, dimension-A: 
continuum operators at point x which are invari- 
ant under gauge transformations and all symme- 
tries of the anisotropic lattice. There are no such 
operators for dimensions less than four. There 

(4) 

are only two dimension-four operators: Q\ = 
g^TrE^ and Q^'*^ g'^TrB'^, where E and B are 
the chromoelectric and chromomagnetic fields, re- 
spectively. There are no dimension-five operators, 
but there are 18 independent dimension-six oper- 
ators, ten of which may be expressed as unimpor- 
tant total derivatives. Thus, we adjust the cou- 
phngs of the W/3[C/] lattice operators so that the 
coefficients of the eight dimension-six operators 
vanish and the coefficients of the two dimension- 
four operators equal each other. If we are inter- 
ested only in on-shell improvement, we are also 
free to use a field redefinition to set some of the 
coefficients of the dimension-six operators to zero. 
Perturbation theory by itself does not reliably 



determine the couplings in the improved action. 
It is known that a judicious combination of per- 
turbation theory and mean field theory works 
much better. Mean field theory is introduced 
by renormalizing the link variables: Uj{x) —^ 
Uj{x)/Us and Ut{x) Ut{x)/ut, where Ut and Us 
are the mean values of the gluon link operators 
for the temporal and spatial links, respectively. 

Carrying out the above procedure, one of_the 
simplest actions we can arrive at is given 



Si[U] = 



tp 



1 



(2) 



where (3 = G/g^, Wc = Y.c Tr(l - [7^), and 
?7sp denotes the spatial plaquettes, C/tp indicates 
the temporal plaquettes, Usr denotes the product 
of link variables about a planar 2x1 spatial rect- 
angular loop, Ustt refers to the short temporal 
rectangles (one temporal spacing, two spatial), 
and Uttr refers to the tall temporal rectangles 
(two temporal spacings, one spatial). The mean 
values Ut and Us are determined by guessing in- 
put values for use in the action, measuring the 
mean links in a simulation, then readjusting the 
input values accordingly. When at is significantly 
smaller than Os, we expect the mean temporal 
link Ut to be very close to unity. For example, in 
Landau-gauge perturbation theory, 1 — (iTrC/t) « 
(ot/as)^. Hence, to simplify matters, we set 
Ut = l. A convenient and gauge-invariant defini- 
tion for Us in terms of the mean spatial plaquette 
is then given by Ua = (|ReTrPss')^^^, where Pss> 
denotes the spatial plaquette. 

One defect of the action Si is that the gluon 
spectrum has spurious high-energy states arising 
from Wttr, which spans two time slices. These 
modes occur at energies of order 2 /at and have 
little effect on simulation results, but they can 
cause problems when applying the variational 
method to extract masses from short-time corre- 
lation functions. These modes may be eliminated 
by relaxing the improvement conditions, requir- 
ing that the coefficients of all dimension-six oper- 



Table 1 

Rcnormalization of the anisotropy. The input 
anisotropy is ^ = at/ as, and ^meas is the measured 
anisotropy (see text). The lattice spacing Us is 
set using the IP — IS" splitting in charmonium 
for SU3 and the string tension for SU2- Slw is 
the improved Liischer-Weisz action[^. 



Action 


a 




Cmcas 


as (fm) 


Slw 


1.9 


1/1 


1.00 


0.397(4) 


Sh 


2.0 


1/2 


0.505(11) 


0.418(6) 




2.0 


1/3 


0.342(10) 


0.372(6) 




2.3 


1/2 


0.507(10) 


0.311(5) 


Si 


2.0 


1/2 


0.505(10) 


0.371(4) 




2.3 


1/2 


0.510(11) 


0.271(3) 


Sll{Us=l) 


3.9 


1/2 


0.40(2) 


0.40(2) 


SiiiSU^) 


0.848 


0.276 


0.286(6) 


0.371(9) 




1.027 


0.351 


0.357(4) 


0.283(6) 




1.114 


0.409 


0.417(6) 


0.249(6) 



ators except Tr[(D(E)^] vanish. One then obtains: 
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This action has 0{a^) errors which are very small, 
being suppressed by a factor of For this ac- 
tion, various values for the parameter Us{P,S,) 



0.745, u.(2.0,i) 



0.772, and 



are Ms(1.7, i) 
w,(2.4, 1) = 0.806 

The rcnormalization of the anisotropy can be 
determined by measuring the static-quark poten- 
tial V{x, y, z) from Wilson loops in different ori- 
entations. For example, Wilson loops in the xt 
and xy hyperplanes have the following asymptotic 
behaviours (/, J positive integers): 



Wxt{Ias,Jat) 
Wxy{Ias, J as) 



J — ^00 



Zxte 



-Ja,[V{Ia,,0,Q) + Vo] ^g-j 



where Z^t, Z^y, and Vq are rcnormalization con- 
stants. From the asymptotic behaviour of Wxt we 
obtain A;^;* a* [V{2as, 0, G)-V{as, 0, 0)] and from 
Wxy we obtain A^^ = a,[y(2a,, 0, 0)-V^(a,, 0, 0)]. 
The renormalized anisotropy can then be defined 
as the ratio ^meas = ^xt/^xy Results for the 
renormalized anisotropy [Q are listed in Table 0. 
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Figure 1. The static-quark potential from Wil- 
son loops in different orientations. Results are 
shown for the action 5// with f3 — 2.0 and 
at /as = 1/3. The different potentials were shifted 
to agree at r = as- 



Note that without tadpole improvement (us — 1), 
the renormalization of the anisotropy is a 20% ef- 
fect; mean- field improvement reduces this renor- 
malization to the 1 — 3% level. Results for the 
static-quark potential[|| are shown in Fig. [|. 

The masses of various low-lying glueballs were 
also calculated The results are shown in Fig. || 
as a function of the lattice spacing as measured in 
terms of the hadronic scale tq. The scalar glue- 
ball mass from the improved action Sjz exhibits 
dramatically reduced cutoff contamination com- 
pared to the Wilson action. Finite-Os errors are 
seen to be small for the tensor and pseudovector 
glueballs, although differences between the 
and T^~^ representations indicate violations of ro- 
tational invariance, especially for large a^. 

These results clearly demonstrate the useful- 
ness of perturbatively-improved gluonic actions 
on anisotropic lattices. For coarse-lattice glue- 
ball studies, anisotropic lattices are crucial. The 
effectiveness of tadpole improvement in reducing 
the renormalization of the anisotropy to the level 
of a few per cent makes using anisotropic lattices 
no more difficult than isotropic ones. This work 
was supported by grants from the NSF, the DOE, 
and the NSERC of Canada. 
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Figure 2. Low- lying glueball masses mg against 
lattice spacing Os in terms of the hadronic scale 
ro « 1/2 fm defined by [r^ dV (r) / dr]r=ro =1.65 
where V{r) is the static quark potential. The x 
denote results for the scalar and tensor glueballs 
using the simple Wilson action||^. Results using 
the improved action Su for the Af^, E~^^ , T^^, 
and representations are indicated by o, □, 

0. and >(!, respectively. The dashed lines are linear 
fits to the Wilson results for the scalar and tensor; 
for the 1^ , the dashed line guides the eye. 
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